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Abstract 

We study the potential of the hi aioi — > 4r signal from the lightest scalar (hi) 
and pseudoscalar (ai) Higgs bosons to cover the parameter space of the Next-to- 
Minimal Supersymmetric Standard Model (NMSSM) at the Large Hadron Collider 
(LHC). We exploit a 2fi -\- 2jets signature from four r's decays (accompanied by 
missing transverse energy), resorting to both Higgs-strahlung (HS), by triggering 
on leptonic decays, and Vector Boson Fusion (VBF), by triggering on two same 
sign non- isolated muons. 
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1 Introduction 



The NMSSM [1-13], which contains a singlet /singhno superfield in addition to the particle 
content of the Minimal Supersymmetric Standard Model (MSSM), has various advantages 
with respect to the MSSM such as the solution of the /i problem and the smaller fine- 
tuninfl [18-20]. The additional singlet field results in one new CP-even and one new CP- 
odd state in the Higgs sector which consists of three CP-even mass eigenstates (/ii,2,3) and 
two CP-odd states (01,2) in the NMSSM. In contrast to the MSSM, it is not certain that 
even at least one Higgs boson can be found at the CERN Large Hadron Collider (LHC) 
within the NMSSM. This problem can be created by the additional singlet states in the 
Higgs sector leading to invisible decay channels of the SM-like Higgs states. Of particular 
relevance for the Higgs search at the LHC in the NMSSM are the decays of the lightest 
CP-even Higgs into lighter CP-odd states, hi aiOi. This channel has been claimed 
to be the only means to establish a no-lose theorem for the NMSSM [21-31], at least 
in parameter regions where the Supersymmetric (SUSY) partners of ordinary Standard 
Model (SM) objects are made suitable heavy. Also notice that the hi state could well be 
below the LEP limit on the SM Higgs mass, of 114 GeV, albeit with weakened couplings 
to ordinary matter. 

In this letter we will focus on NMSSM parameter regions with light ai states (light ai 
scenario) with M^^ < 2mj, where the decay ai — > bb is not possible. In these parameter 
regions with large BR(ai — > r+r^) the scope of hi aiOi decays into jjr^r^ pairs (where 
j represents a jet of either heavy or light flavour and where the r's decay leptonically) 
has been found to be rather questionable [32]. Hence, here we investigate the scope of 
the 4r channel, wherein two r's are searched for in their muonic decays, while the other 
two are selected via their hadronic ones. In doing so, we exploit both the HS and VBF 
production channels. 

The plan of the paper is as follows. In section 2 the results of a detailed parameter 
scan are discussed identifying the parameter regions with light ai and possibly light hi 
with M/i^ < 114 GeV. In section 3 the hi production and decay processes are simulated 
with the PYTHIA event generator [33] and the prospects for the determination of the 
NMSSM Higgs bosons in the channel h At are discussed. We then conclude 

in section 4. 

2 The low-energy NMSSM parameter space for the 
hght ai scenario 

Here we investigate the NMSSM parameter space for Ma-^ < 2mb, with particular interest 
in the cases where the aforementioned hi —>■ aiai 4r decays may be visible at the LHC 
from HS and/or VBF production processes. The NMSSM Higgs at the Electro- Weak 
(EW) scale is uniquely defined by fourteen parameters: the ratio of the doublet Higgs 

^Alternative formulations of the MSSM with an additional singlet/singlino superfield — known as 
the Minimal Non-minimal Supersymmetric Standard Model (MNSSM) and new Minimally-extended Su- 
persymmetric Standard Model or nearly-Minimal Supersymmetric Standard Model (nMSSM) — also 
exist [14-17]. 
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vacuum expectation values (VEVs) tan/?, the trilinear couplings in the superpotential 
A and k, the corresponding soft SUSY breaking parameters A\ and A^, the effective 
parameter fi = X{S) (where (S) denotes the VEV of the singlet Higgs field), the gaugino 
mass parameters Mi, M2 and M3, the squark and slepton trilinear couplings At, Ab and 
At, and the squark and slepton mass parameters Mj^ and Mf^. In the following we 
will establish the NMSSM parameter regions, defined in terms of the above inputs, that 
survive present theoretical and experimental constraints. 

2.1 Full NMSSM parameter scan 

First we perform a 'wide' scan over the full NMSSM parameter space, where the numerical 
values over which the parameters have been varied are: 

-1000 GeV <A^< 100 GeV, - 5 TeV < Ax < 5 TeV, 100 GeV < fi< 1000 GeV, 

10"^ <X,K< 0.7, 1.5 < tan/5 < 50 (1) 
while the remaining parameters (entering the Higgs sector at loop-level) were fixed at 

M1/M2/M3 = 150/300/1000 GeV, At = Ab = A^ = 2.5 TeV, Mj^ =Mj^ = l TeV, 

(2) 

resulting in a heavy spectrum of the SUSY partners oi O = 1 TeV (except for the U(l) 
and SU(2) gauginos). This scan (as well as all the others in the remainder of this letter) 
has been performed by using the NMSSMTools package [34-36], which calculates NMSSM 
spectra (masses, couplings and decay rates) and takes into account experimental inputs 
including LEP limits, -B-physics bounds as well as constraints on the relic density of 
cold dark matter (DM) stemming from the results obtained by the Wilkinson Microwave 
Anisotropy Probe (WMAP) [37]. 

In Fig. [T] we present the results of this scan in form of scatter plots, corresponding to 
10^ generated points, where about 2.7 x 10^ points survive all constraints but the DM 
constraints (red points) and about 7.6 x 10^ points survive all constraints (black and green 
points). The bands between 115 GeV < Mh^ < 130 GeV appear because in a majority 
of the generated points the hi has SM-like properties where the lower bound of the band 
is given by the LEP Higgs limit and the upper bound by the theoretical upper limit of 
the hi mass in the NMSSM. For the points with M^-^ ^115 GeV the hi couplings to SM 
gauge bosons are suppressed, hence the LEP Higgs limit is not applicable. Though only 
17 points with Ma^ < 10 GeV survive (green points), one can see from Fig. [^a) that 
small |Ak|'s are preferred while Fig. [H^b) indicates the preference for large positive A\. 

2.2 Scan for narrowed ^4^ range 

The results of Fig. [1] (especially the preference for small A^'s) motivated us to 'narrow' 
the range of the parameters, by scanning it over the intervals 

- 20 GeV < A, < 25 GeV, - 2 TeV < < 4 TeV, 100 GeV < < 300 GeV, (3) 

and the rest of the parameters as in Eq. ([T]). Fig.O based on 10'' generated points resulting 
in 2.9 X 10^ points surviving all but the DM constraints (red points) and 4.0 x 10^ points 
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Figure 1: Results of the 'wide' scan over the full NMSSM parameter space, Eq. ([T]), 
mapped onto the planes: (a) [A^,Mh^], (b) [Ax,Mh^], (c) [k, M/^J, (d) [A,M,jJ, (e) 
[Mq^,M/ij], (f) [Maj,m^o]. Colour code: red - all constraints are satisfied but rehc 
density (above WMAP constraint: VlK^ > 0.11); black - all constraints are satisfied, 
Mai > 10 GeV; green - all constraints are satisfied, Mq^ < 10 GeV. 

surviving all constraints (black and green points), shows that this is precisely the region 
where a large portion of the NMSSM parameter points with M^^ < 10 GeV (here 3.3 x 10'^ 
green points) are found. Now we can also see certain correlations and other interesting 
features onsetting in the Ma^ < 10 GeV region: 

1. values of Ax > are preferred, see Fig. [2^b); 

2. for points with low Mh^ small values of (Fig. EJ^a)) as well as small values of k 
(Fig. [2](c)) are preferred; 

3. in the Mq^ < 10 GeV region of our interest, Mh^ as low as 20 GeV can appear 
(Fig.EKe)). 



2.3 Final scan for the light ai scenario 

We have then performed a 'final' scan over the NMSSM parameter space by requiring 
Mai < 10 GeV additionally to Eq. ([H]). At this stage, we recast the NMSSM parameter 
space in terms of the quantities entering hi HS and VBF production and hi aiai 4r 
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Figure 2: Results of the scan over the 'narrowed' NMSSM parameter space, i.e., analogous 
to Fig. □ but for -20 GeV < A« < 25 GeV, -2 TeV < Ax < A TeV, 100 GeV < /x < 
300 GeV. The individual plots and the colour code are the same as in Fig. [TJ 

decays. The results of this scan (10^ generated points, resulting in 3.9 x 10^ points with 
Mai < 10 GeV surviving all constraints) are shown in Fig. [3l where the colours were chosen 
to indicate the measure of decoupling of the lightest CP-even Higgs boson, hi, from the SM 
limit (denoted simply by H). This measure of decoupling was defined through the ratio of 
the coupling strength (squared) of the ZZhi vertex in the NMSSM relative to the SM case 

(in fact, this is the same for couplings to gauge bosons): Rzzh = (s'zzhf ^/S'Izh) )• 
One should notice that both HS, pp — > Vhi, and VBF, pp — * jjV*V* — > jjhi, rates 
{V = Z, W"^) are directly proportional to Rzzh and are suppressed in the non-decoupling 
regime whenever Rzzh is essentially smaller then unity. 

From Figs. [3] and m one can see the following important features of the Ma^ < 2mh ~ 10 
GeV scenario: 

1. The lighter the Higgs hi the more significant should be the NMSSM deviations 
from the SM case, e.g., for any M^^ < 50 GeV any Rzzh is limited to be < 0.5, as 
dictated by LEP constraints [38, 39] (this correlation is illustrated in a more clear 
way in Fig. 111(a), presenting the Rzzh versus Mh^ plane, which exhibits the typical 
pattern of the LEP Higgs exclusion curve [39]). 

2. In the Mu^ < 40 GeV region Ax is always positive (Fig. [^b)), A < 0.45 (Fig. [3](d)) 
while K < 0.1 (Figs. [3](c), m^b)), which corresponds to an approximate Peccei-Quinn 
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Figure 3: Results of the 'final' NMSSM parameter scan, i.e., with M^^ < 10 GeV and 
with Eq. (IHl) enforced. The black, blue and green colours indicate the cases Rzzh < 0.1, 
0.1 < Rzzh < 0.5 and Rzzh > 0.5, respectively (where Rzzh is defined in the text). The 
individual plots are the same as in Fig. [H 

(PQ) symmetry (/? ^ and —>■ 0) [13,40]. However, for M^^ > Mz the whole 
range < k < Kmax ~ 0.7 is populated. Hence an approximate PQ symmetry is 
not necessary for the M^^ < 2mh ^10 GeV and M^^ > Mz scenario while some 
fine-tuning ensuring l^^l <^ \Ax\ is present. 

3. One should further notice the correlation between the singlet nature of the hi and the 
singlino component of the lightest neutralino as well as the correlation between their 
masses, which are visualised in Fig. |l](c) and (d), respectively. These correlations 
take place to satisfy WMAP constraints on relic density. In fact, in the lower- 
left corner Fig. Hl^d) we can clearly observe XiXi ~^ hi the annihilation corridor 
for M/jj < 80 GeV, along which 2m^o ~ M^-^ that allows one to lower down the 
neutralino relic density to a level consistent with experimental constraints. In this 
corridor Xi and hi have significant singlino and singlet components, respectively. 
At the very bottom of Fig. IH^d) we can further see another DM motivated region 
with very light (< 5 GeV) neutralino where 2m^o ~ Ma^ provides the XiXi ~^ (^i 
annihilation corridor. This region is also represented in Fig. [3]^f) as a border-line 
at the left edge of the parameter space. These two representative DM motivated 
regions appear because the NMSSM structure requires hi to be a singlet and Xi to 
be a singlino for Ma-^ < 10 GeV and low M/^^. 
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Figure 4: Results of the 'final' NMSSM parameter scan, i.e., with M^^ < 10 GeV and with 
Eq. ([3]) enforced: (a) Rzzh (see text) versus M/j^, (b) k versus A^, (c) Rzzh versus the 
singlino component yl^ of the lightest neutralino x?, (d) m^o versus M/^^, (e) BR(/ii — > 
OiOi) versus M^^, (f) BR(ai t^t^) versus M^^. The colour coding is the same as in 
Fig. El 

4. The branching ratios BR(/ii ^ aiOi) (Fig. H^e)) and BR(ai r+r") (Fig. H^f)) 
are large for a large fraction of the parameter points with < 2mf, ~ 10 GeV 
which is crucial for the hi — > OiOi — > 4r mode analysed in the next section. 

3 Phenomenology of the hght ai scenario 

As a final step of our analysis, we have combined the production rates of VBF and HS with 
selection efficiencies evaluated by generating these processes within the PYTHIA Monte 
Carlo (MC) generator. The latter have been estimated in presence of experimental-like 
cuts, after parton showering and hadronisation and with underlying events turned on. The 
jets, in particular in the VBF analysis, were found with the PYTHIA routine PYCELL 
using the approximate CMS calorimeter tower granularity and cone size 0.5. No detector 
smearing of the jet energy is applied, however, we do not expect a significant reduction 
of the event selection efficiency once this is enforcecl^. 

The hi —>■ aitti t~^t~t~^t~ mode from VBF hi production has been considered for 



^The selection efficiency is related to a cut on jet transverse energy Et and since we cut right on 
the peak of the tagging jet Et distribution we do not expect jet energy smearing to affect this efficiency 
substantially. 
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the /i^/i^Tj^^Tj^^ final state containing two same sign muons and two same sign r jets (r 
jet means the r lepton decaying hadronically, r — > hadrons + u) . The r leptons from the 
light ai decay are approximately coUinear. Fig. [5](a) shows the separation in the {rj, cj)) 
space between the two r leptons and between the muon and the r jet in the decay chain 
ai TT —>■ fiuu + Tjetz^ for the benchmark point P2, M/^^ = 120.2 GeV and = 9.1 
GeV, proposed in [18,20]. The non-isolated, di-muon High Level (HL) trigger is needed 




Figure 5: (a) Separation in (77, 0) space between the two r leptons and between the 
muon and the r jet from the ai tt fiuu + TjetZ^ decay chain in the VBF hi 
OiOi — > T^T~T~^T~ fi^^^T^^^T^^ channel for the benchmark point P2: M/i^ = 120.2 GeV, 

Ma-^ = 9.1 GeV. (b) Distributions of p'^, and p^'^'^^ for the point P2. 



to select the signal events. The standard CMS di-muon trigger with relaxed isolation has 
a di-muon threshold 10 GeV on both muons for C = 2x10^^ cm~^s~^ [41]. The muons 
from the signal events are very soft as one can see in Fig. [5](b), thus lower thresholds are 
needed. For example with the 7 GeV threshold the efficiency is increased approximately 
by a factor of two. However the QCD background rate is also increased by approximately 
a factor of two [42] which is not acceptable. In order to cope with the rate, the same 
sign relaxed di-muon trigger was introduced recently in the CMS trigger table [43]. The 
rate of di-muons from bb production is reduced by a factor of four when asking for two 
muons of the same sign at the threshold 5 GeV. The off-line selection strategy requires 
the presence of the two same sign, non-isolated muons with one track within a cone 0.6 
around each muon direction, thus selecting the one prong r decays. The full list of off-line 
selections is: 

• two same sign muons with px > 7 GeV and |?7| < 2.1 with one track of p-r > 2 GeV 
in a cone 0.6 around each muon. The muon and the track should have opposite 
charge; 

• two r jets with px > 10 GeV, \ri\ < 2.1; 
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pp ^ Wh^^ N a^a^ /v4t rates in fb pp jj Aj, jj a^a^ jj 4% rates in fb 




M,,(GeV) M,,(GeV) 

Figure 6: Cross sections (including all relevant branching ratios) for HS (a) and VBF (b) 
after the selection cuts described in the text. The population of points used correspond to 
that of the 'final scan' described previously. Black diamonds correspond to the benchmark 
P2 (right) and P3 (left) from Refs. [18,20]. 



• two jets with pj- > 30 GeV, \ri\ < 4.5. 

The hi —>■ aitti —>■ t^t^t^t~ mode from HS production with leptonic decays of the W 
bosons can give a very clean, almost background free signal. The lepton coming from the 
W decay is used for the triggering of the event using the CMS single isolated lepton trigger 
with the thresholds 19 GeV and 26 GeV, respectively, for the muon and electron [42]. The 
final state £ + /xTjet/^Tjet (£ = e or /i, where the lepton comes from W decay) has been 
considered and r jets with one track have been selected. For each r-jet candidate there 
must be a muon in a cone 0.6 around the track. Unlike the VBF case, the muons are not 
required to have the same sign. The full list of selections is: 

• trigger selection: isolated muon or electron with px greater than 19 or 26 GeV, 
respectively, and \rj\ < 2.5; 

• two 1-prong T jets with px > 10 GeV, |?7| < 2.1. For each r jet one muon of px > 7 
GeV within a cone 0.6 around the r track should exist and have a charge opposite 
to the track charge; 

• events with extra jets in addition to the two r jets are rejected. 

The results in Fig. [6]show that, after our final scan, the population of parameter points 
is such that in both channels the highest cross sections are found for M^^ > Mz, although 
in the case of VBF also lower hi masses can yield sizable rates. Independently of M/j^, 
the Oi mass enables sizable event rates anywhere above 2m^, but particularly just above 
the threshold. At high luminosity, 100 fb~^, the highest rates would correspond to 1000 
events per year for HS and 8000 for VBF. The black diamond symbols in Fig. [6] denote the 
two NMSSM benchmarks points P2 (right point) and P3 (left point), defined in [18,20]. 
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As it can be appreciated, they correspond to event rates that are mid range amongst all 
those explored, hence not particularly biased towards a far too favourable NMSSM setup, 
yet susceptible to experimental discovery. 

4 Conclusions 

We have analysed the CP-even Higgs boson hi in the NMSSM decaying into aiOi pairs in 
turn yielding four r leptons, produced at the LHC in HS and VBF production and searched 
for through their semi-leptonic/hadronic decays into muons and jets. We have found that 
there is significant potential, especially via the VBF production channel, to discover a 
Higgs boson at the LHC in NMSSM scenarios with M^^ < 2m5 which is an important step 
in establishing a no-lose theorem for the NMSSM. In order to achieve this the same sign di- 
muon trigger with lower threshold [43] is crucial. We have restricted ourselves to the case 
< 2m(,, where the hi decay fraction into r's in enhanced (otherwise Oi — > bb decays 
are dominant). A scan of the low-energy NMSSM parameter space without assuming 
unification at the high scale has shown that the hi state can be very light, indeed at times 
lighter than the ai. This configuration docs not take place in the constrained NMSSM 
(cNMSSM) of Ref. [18, 19]. However, we are currently investigating whether such light 
hi masses can be found in the NMSSM with non-universal boundary conditions at the 
unification scale [44]. With reference to the NMSSM benchmarks points defined in [18,20], 
we note that those relevant to our 4t channel are the P2 and P3 benchmarks, which yield 
event rates in the mid range amongst those explored here. Our summary is preliminary, 
as only signal processes have been considered and only in presence of MC simulations, 
with no backgrounds and full detector performance enabled. The latter clearly ought 
to be investigated before drawing any firm conclusions and this is currently being done. 
However, we would like to conclude that the 2/x -|- 2jets signature from hi ciicii ^ 4t 
decays with the hi produced in HS and VBF processes does produce sizable event rates as 
high as 1000 events per year for HS and 8000 for VBF so that we expect that a significant 
part of the viable NMSSM parameter space will be covered by using the signature and 
selection cuts suggested in our study. 
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